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Abstract. We observed the quadruple gravitationally 
lensed image of MG J0414+0534 on 23 November 1997 
with a global VLBI array at 8.4 GHz. We report wide-field 
imaging results of its four components at submilliarcsec- 
ond resolution, displaying complex core-like and jet-like 
extended structures. A simple model combining a singu- 
lar isothermal ellipsoid to represent the main lens galaxy, 
external shear, and a singular isothermal sphere to rep- 
resent an additional, nearby object accounts well for the 
core positions and flux densities of the VLBI images. This 
model predicts delays between the different lensed images 
of several weeks. 

Key words: Techniques: interferometric — Galaxies: jets 
— Galaxies: quasars: individual: MG J0414+0534 — Cos- 
mology: gravitational lensing 



1. Introduction 

The gravitationally lensed image of QSO MG J0414+0534 
(z = 2.639 ±0.002 (Lawrence et al. p.995[ ), V = 21.22, also 
PKS 0411+05, 4C 05.19, and DA 131) with four compo- 
nents separated by up to 2 arcseconds is one of the best 
examples of morphology produc ed by an elliptical gravi- 
tational potential (Hewitt et al. 1992 , Katz et al. 1997|) . 
The le ns ap pears to be an elliptical galaxy (Schechter & 
Moore |l993J), a t a redshift z = 0.9584 ± 0.0002 (Tonry & 
Kocliaiiek~ Tl999| ). 

Different studies in radio using the VLA have pro- 
vided important information about the structure of the 
gravitational lens (Moore & Hewitt ( |l997j ) and refer- 
ences therein). Hubble Space Telescope (HST) observa- 
tions show an optical arc joining the three brightest com- 
ponents (Falco et al. 1997). VLBI images of this radio 
source show some extended structure in the sub-images 



(Patnaik & Porcas 1996 Porcas 1998) 
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We have undertaken a project of multi-epoch observa- 
tions of the quadruple images of each of the gravitationally 
lensed objects MG J0414+0534 and B1422+231, to com- 
pare possible structural changes in the lensed images and 
measure possible shifts in the relative positions of the com- 
ponents. These shifts would be caused by intrinsic changes 
in the lensed object or by an alteration of the alignment 
in the system observer-lens-lensed object due to relative 
pr oper m otion of the latter two elements (see Kochanek et 
al. 1996 ). For non-magnified extragalactic objects and pe- 
culiar velocities of, say, 100 to 1000 km s^ 1 , such relative 
proper motions would not be detectable astrometrically 
over decades. Because of their high degree of symmetry, 
quadruple gravitationally lensed systems can magnify in- 
trinsic velocities by factors of 10-100 while minimizing 
systematic effects for measurements of relative motions 
(Kochanek et al. |1996| ). 

We present here the mapping and modeling results 
from the first epoch of observations of MG J0414+0534. 
We describe our VLBI observations in Sect. 0, report our 
imaging results of MG J0414+0534 in Sect. |T discuss our 
modeling results in Sect. Hand give a summary in Sect, pi 

2. Observations and data reduction 

We observed MGJ0414+0534 at 8.4 GHz on 23 Novem- 
ber 1997 with a global VLBI array for over 12 hr. We 
used a basic integration time of 0.4 s. The antennas used 
were (name, diameter, location): Medicina (32m, Italy), 
Noto (32m, Italy), Effelsberg (100m, Germany), DSS63 
(70m, Spain), the Very Long Baseline Array (VLB A) 
(10 antennas of 25m across the US), the phased- VLA 
(130m-equivalent, NM, US), and DSS14 (70m, CA, US). 
Data were recorded with right circular polarization us- 
ing mode 128-2-2, achieving a bandwidth of 32 MHz, and 
correlated at the National Radio Astronomy Observatory 
(NRAO) VLBA Correlator (Socorro, NM, US). The ob- 
serving mode put the correlator close to most demanding 
performance, taking advantage of its latest technical im- 
provements (actually, achieving the smallest sample inter- 
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val -widest channel bandwidth- with the highest spectral 
resolution -largest delay range- in the used mode). The 
antennas pointed to the center of the field of the four im- 
ages of the radio source (all within the individual antenna 
beam at the same time) and only one correlation pass was 
needed to determine the interferometric visibilities jointly 
for the four components. Due to operational reasons, the 
strongest of the four components (Al, see below) was cho- 
sen as phase center at correlation time. 

The observing duty cycle consisted of 2 min on the 
calibrator source TXS 0357+057 and 5.5 min on the tar- 
get source MGJ0414+0534. We allocated some calibra- 
tors and fringe-finders (BL 1803+784, QSO 1928+738, 
BL 2007+777, 3C 84 and 3C 345) at the beginning and the 
end of the observations. 

We analyzed the data using the Astronomical Image 
Processing System (aips, version 15oct99). We applied 
a priori amplitude calibrations determined from the mea- 
sured system temperatures and available antenna gains. 
We aligned the phase slopes across each of the 16 MHz 
bands using the information from each telescope, when 
available. In the absence of this information, we applied 
manual phase calibration deduced from observations of 
TXS 0357+057, 3C84, and 3C345. We then performed 
standard fringe-fitting on the calibrators and the target 
source. We made a preliminary map of MG J0414+0534 
by down- weighting (tapering) the longest baselines to gain 
sensitivity to extended surface brightness. We used this 
initial image to perform a second fringe-fit data correc- 
tion in MPS. 

Since the total angular size of MGJ0414+0534 is 
2 arcsecond s, wid e-field imaging techniques (see, e.g., Gar- 
rett et al. 1999 ) are needed. At a given frequency, the 
field of view of a VLBI map is determined by the av- 
eraging time of the data and th e bandwidth of t he ob - 
servation (see, e.g., Cotton 199S , Bridle & Schwab 1999| ). 
Over-averaging in frequency and time leads to bandwidth- 
and time-smearing, respectively. For our data, we aver- 
aged every four 0.25 MHz channels at the two 16 MHz fre- 
quency bands to obtain 32 frequency channels, each of 
1 MHz width. Later on, to speed up the mapping pro- 
cess, we time-averaged in 4 s intervals (shorter averaging 
times, 1.2 and 2.4 s, did not produce significant varia- 
tions). The adopted 4 s and 1 MHz averagings avoid both 
types of smearing, allowing a size of the data set just small 
enough to be manageable with our computing resources. 
This wide-field mapping technique allows to image fields of 
arcsecond-size with milliarcsecond resolution. However, a 
single mapping field results of unmanageable size. As a so- 
lution, the AIPS task IMAGR provides a reliable wide-field 
mapping methodology based on a multi-windowing ap- 
proach. Our imaging strategy emphasized tapering at the 
beginning of the mapping process to help recovering most 
of the extended flux density in the sub-images. In further 
iterations we lifted the tapering to achieve the full reso- 
lution provided by the global VLBI array. Using a hybrid 
weighting (between uniform and natural, robust=— 1 in 



IMAGR task) at full resolution, the synthesized beam ob- 
tained was of 2.55xl.l3mas (position angle (P.A.) 14.6°). 

3. Imaging results 

We show the final AIPS maps in Fig. [y. At the center of 
Fig. [l| we plot a wide-field image of the source, restored 
with a circular Gaussian beam 25 mas in diameter to il- 
lustrate the location of the four components. At this reso- 
lution, MG J0414+0534 comprises two bright, close com- 
ponents (Al and A2) separated by 0''47 and two weaker 
components B and C separated 2" 12 from each other, and 
2!'02 and 2''08 from Al, respectively. We show the map 
dimensions in the panel corresponding to Al (bottom, 
left), preserving the same map scale in the other three 
sub-images. Each of these sub-images has been restored 
with the appropriate synthesized beam, which is about 
2.55 x 1.13 mas, P.A. 14.6°. In Fig. g, we show more de- 
tailed maps of the brightest features of the four sub-images 
at full resolution. We present in Table [j] their peaks of 
brightness / pea k and the total flux density values (Stat) for 
the whole sub-image and for core-like and jet-like struc- 
tures, obtained computing the flux densities of the CLEAN 
components in each of those regions. VLA flux density 
ratios among sub-images given by Katz & Hewitt (1993) 
and Katz et al. ( 1997| ) are similar to ours. Our results show 
Al and A2 somewhat brighter, and B somewhat weaker 
as their results. It is difficult to compare both data sets, 
since some flux density present at the VLA may not be 
detected in the VLBI images. 

Table 1. Total flux density values for each of the sub- 
images and their sub-components in Fig. |l|. 
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Relative position of the peaks of brightness on the maps. 

The image pixel size is 0.16 mas. 

Total flux densities in the fields shown in Fig. |y. 

Total flux densities for the complex core-like region shown 

in Fig. |. 

Total flux densities for the jet-like structures, labelled as 

jet and c(ounter)-jet. 



The gravitational lensing geometry produces four im- 
ages of the same object, deformed and with different mag- 
nification factors. In addition, some lensed, unresolved 
components might also be fainter than the sensitivity level 
of our maps, causing some confusion in the feature identifi- 
cation among different sub-images. In our maps, we do not 
observe isolated, compact components (like in the quadru- 
ple lensed object B1422+231) but, rather, very complex 
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Fig. 1. Composition of our images of the four gravitationally lensed images, labelled as Al (bottom, left), A2 (top, left), 
B (top, right), and C (bottom, right) of MG J0414+0534 at 8.4 GHz with a global VLBI wide-field image (central panel) 
of the radio source. The axes are relative right ascension and declination, in milliarcseconds. The wide-field image rep- 
resents a combination of the four mapped windows (or sub-images) in a wider grid convolving their CLEAN-componcnts 
with a circular beam of 25 mas in diameter. It displays contours of 1 mJy/beamx (—1, 1, 2, 4, • • •). The sub-images are 
convolved with their synthesized beam of 2.55xl.l3mas (P.A. 14.6°). The contours are 0.5mJy/beamx(— 1, 1, 2, • ■ •). 
The root-mean-square (rms) noise levels in the sub-images were 0.25, 0.21, 0.16, and 0.16mJy/beam, for Al, A2, B, 
and C, respectively. The values of the peaks of brightness of each of the sub-images are given in Table W. 



structures, with core-like and jet-like regions that show aligned with respect to the core-like region. The fainter 



different morphologies for each lensed image. 

We discuss in turn the features of each lensed im- 
age. Image Al. It is the most complex and the brightest. 
There are two parallel jet-like components, totally mis- 



one is 20 mas away, P.A. ~ —15°. The brighter one, very 
stretched along P.A. ~ 30° begins at about 10 mas south- 
west of the core-like region, and extends over more than 
50 mas. The brightest feature (core-like) has two compo- 
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Fig. 2. Sub-image cores con- 
volved each with the appro- 
priate beam ~2.55xl.l3mas, 
P. A. 14.6°. The axes are relative 
right ascension and declination, 
in milliarcseconds. Contours are 
0.6mJy/beamx(-l,l, y/2,2,---). a 
and b in the maps indicate the two 
brightest features in each sub-image 



nents along P.A. ~ —60° and a resolved structure. Im- 
age A2. Two extended regions not quite co- aligned (in 
the same direction) are present at about 30 mas to the 
NE and 80 mas to the SW of the central (core- like) area. 
The core-like area (Fig. g, top, left) is symmetric, with 
two components, elongated over more than 15 mas at P.A. 
~ 20°. Image B. It presents two extended jet-like regions, 
fairly aligned with its central area, at about 70 mas to 
the NW and SE. The core-like area itself is double and 
elongated more than 15 mas at P.A. ~ —70°. Image C. 
It presents a double N-S structure, the northern feature 
being the brightest, and extended emission to the north of 
this double structure. The core-like region presents some 
extended emission to the South. 

A comparison of our maps with previously published 
images of MG J0414+0 534 sh ows that the EVN images at 
5 GHz shown in Porcas (1998 ) are similar to our maps, but 
with a wider beam and lower dynamic range. The same 
holds for the 1.7 GHz EVN images from Patn aik fc Porcas 
( |1996| ). The optical image from Falco et al. ( |1997| ) shows 
emission extending over an arc joining Al, A2 and B. 



4. Lens models 

Quadruple gravitational lens systems can be produced by 
an elliptical lensing potential. Many attempts have been 
made to model MG J0414 +0534 (e.g. Kochan ek|l99l| , He- 
witt et al. |1992| Witt et al. |l995| ). Falco et al. ( |l997[ ) could 



not find models that were simultaneously consistent with 
the brightest features of the optical images and the lower- 
resolution radio data from the VLA (e.g., Katz & Hewitt 



1993 ). The image flux densities at radio wavelengths have 



shown little variation, and have not yielded time delays 
(Moore & Hewitt 1997). We show below that the addition 
to the lens model of a previously known but unsuspected 
lens galaxy near the images suffices to account for the 
properties of the images. This new model also predicts 
time delays between the different lensed sub-images. 

We fitted a simple lens model to our VLBI data. Be- 
cause the lens galaxy G is not detected in these data, we 
used the deep HST NICMOS H-band (A = 1.6 ^m) im- 
age from the CfA/Arizona Space Telescope Lens Survey 
(CASTLES; e.g., Falco et al. |l999j ) to estimate the coordi- 
nates of the center of brightness and the structural param- 
eters of the lens galaxy G (see Fig. ||) . We registered the 
VLBI and point-like infrared images by assuming that the 
infrared positions of Al, A2, B and C correspond to the 
peaks of brightness of each of the four VLBI sub-images. 
The transformation that maps the infrared coordinates to 
those of a or b is a scale change (to take into account 
a chip position and wavelength-dependent uncertainty of 
~ 10 mas in the NICMOS/NIC2 scale) with a rotation (to 
take into account the uncertainty of ~ 0? 1 in the orienta- 
tion of HST). We found the transformation to the feature 
a (b) requires a rotation of 0?28 (0?16) and a scaling by 
0.992 (0.991), with rms residuals of 11 (5) mas. The radio 
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Fig. 3. CASTLES HST/NICMOS image of 
MGJ0414+0534 at H-band (A = 1.6 /mi), after de- 
convolution of the point-spread function. 

and infrared emission do not emanate necessarily from the 
same region of the source. Therefore, we adopt the bright- 
est feature a as the 8.4 GHz counterpart of the quasar, and 
assign to it positional uncertainties of 11 mas, to reflect the 
uncertainty in the registration. In this paper, we ignore the 
lensed jet structure because of the difficulty in establishing 
the correspondence between lensed components along the 
jets. After ongoing software developments are completed, 
this aspect will be addressed in a forthcoming paper. 

Taking as origin of coordinates the feature a on image 
Al, we use as constraints for the lens models the 6 relative 
coordinates of feature a (Table [l]) (assuming a conserva- 
tive 0.5 mas uncertainty in the positions), the coordinates 
of G, xq = — 1072±llmas, yo = +665±llmas, and the 3 
peak of brightness ratios from the radio sub-images (with 
10% uncertainties due to the possibility of microlensing) . 
We adjusted lens models with a small number of parame- 
ters. Our initial model for the mass density distribution is 
a singular isothermal ellipsoid (SIE), with 5 parameters: 
the Einstein ring radius b, the ellipticity e, sky position an- 
gle 6 e , and coordinates x and y. SIE models yield values 
that are consistent with the masses of other lens galaxies 
(e.g., Lehar et al. 2000D and that reflect the properties of 
elliptical galaxies (e.g., Kochanek et al. 200C). As a rule, 



single mass distributions require external lensing shear to 
account for the properties of quadruple lens systems (Kee- 
ton et al. 1997). Therefore, we added external shear with 



amplitude 7 and position angle Q 1 to our model. The total 
number of parameters of our SIE+7 model is 7. 

We adjusted all 7 parameters of the SIE+7 model us- 
ing the sum of weighted rms residuals for the peak of 
brightness ratios, the coordinates of G and those of the 
a feature images as the figure of merit, x 2 . The best- 
fitting parameters of this initial model yielded \ 2 ~ 30 



principally because the lens fell ~ 60 mas away from 
the observed galaxy. Therefore, we included object X 
(Schechter & Moore |1993| ) at x x = -1457 ± 13 mas, 
y x = +2122 ± 13 mas (derived from the CASTLES data) 
as a possible additional component of the lens. Because 
X is about 10 times fainter than G, its structure can- 
not be distinguished from a circular brightness distribu- 
tion. Thus, we simply assumed X is a singular isothermal 
sphere (SIS) and added its 3 parameters, Xx, yx an d bx 
to the model, where xx and yx are constrained by the 
CASTLES data. We adjusted the resulting 10 parame- 
ters of the SIE+7+SIS model, which therefore has 3 de- 
grees of freedom (DOF). We thus found values that ac- 
count nearly exactly for all the positions of the images 
of feature a as well as those of G and X, and also for 
the ratios of the peaks of brightness at the core-like re- 
gions, with a final merit per DOF x 2 /3 = 1.7 (arising 
solely from the constraints in the peak of brightness ra- 
tios). The parameters of this model are shown in Table 0. 
The model predicts that B is the leading image when the 
source varies: the time delays in a fiat cosmology (fio = L 
H = eSkms^Mpc" 1 ) are At A1B = 15.7 ± 1.3 days, 
AiA2S = 16.0 ±1.4 days and Atcs = 66 ± 5 days, where 
the error bars are an estimation from the variation of the 
predicted time delay as the lens model parameters are var- 
ied within the A% 2 < 1 region. 

Our simple lens model reproduces successfully the 
data. It is remarkable that: a) the position angle of the SIE 
model is consistent with that of the light of G (—96° ± 4); 
b) the ratio of the two Einstein ring radii, which scales 
as the square root of the ratio of luminosities, bx/bc = 
0.17 ± 0.02, is consistent with the brightness ratio of X 
and G, whose magnitudes differ by 2.6 ± 0.5 mag; and c) 
the amplitude of the shear is also consistent with that ex- 



pected from cosmic shear (Bar-kana 1996), and with the 
number of galaxies seen within a radius of ~ 1 arcmin of 
the lens system (McLeod et al. [2000). 



5. Summary 

We present detailed state-of-the-art global VLBI images of 
the gravitational lens system MG J0414+0534. The indi- 
vidual images, separated by distances up to 2 arcseconds, 
exhibit radio structures extending up to 100 mas. The use 
of wide- field mapping techniques in very sensitive observa- 
tions provides high-quality images of extended structures 
in MG J0414+0534, showing features that extend more 
than 1500 beam widths across in the E-W direction, and 
more than 750 N-S. 

We reproduce successfully the relative positions and 
peak of brightness ratios of the radio cores with a lens 
model consisting of a singular isothermal ellipsoid for the 
galaxy G, an external shear to it, and a secondary po- 
tential given by a singular isothermal sphere associated 
with object X. The model predicts that B is the leading 
image and C lags all the other images with delays for all 
image pairs of several weeks. It is important to notice the 
alignment of the optical images of the lens galaxy G and 
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Table 2. Lens model. 
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-83±6 0.096±0.012 56±4 



Positional offsets are relative to the peak of brightness in the Al image. 



Angles from North through East. 
Note: Parameter uncertainties were estimated from the range over which A\ 2 < 1 as each parameter was varied and the 
remaining parameters were optimized. 



the model's elliptical potential (SIE). The model's shear 
is also consistent in amplitude with the expected cosmic 
shear and from the observed environment (within a radius 
of ~ 1 arcmin) of the lens system (McLeod et al. 2000J ) . 
Our model also provides consistency between the lumi- 
nosity ratio for X and G and the ratio of their Einstein 
radii. We plan to refine and study further our lens model 
by adding the constraints provided by the lensed jets, us- 
ing a new software package for lens modeling developed 
by Keeton et al. (in preparation) that avoids the problem 
of matching components along these jets, in the different 
lensed images. 
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